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■ Abstract 
T— I ■ 

I Working within the quasi-metric framework, we examine the gravitational field 

^1 of an isolated, flat disk of spinning dust assuming a weak gravitational field so that 

D I relativistic approximations can be used. Contrary to Newtonian theory, these weak 

Q [ field approximations involve a background 3-geometry being that of the 3-sphere 

^ ! rather than Euclidean 3-space. 

We solve Laplace's equation approximately for small distances from the disk and 
CJ , find two discrete spectra of solutions involving Legendre polynomials and Legendre 

^ . functions, respectively. These spectra correspond to two different mode surface 

bX), density spectra. The physical surface density of the disk can be expanded in terms 

of the surface density mode spectrum involving Legendre polynomials by means of 

CN ■ 

^ . a finite Legendre transform. Moreover, boundary conditions single out the expo- 

. nential disk as a "preferred" physical surface density profile. 

I ■ On the other hand, the surface density modes involving Legendre functions may 

^ ■ also be summed up to yield a second, "induced" surface density (representing col- 

■ lective properties of the disk) that turns out to be singular and thus unphysical. 
^ ■ However, it can be used to construct a non-singular "induced associated" surface 

■ density and a corresponding "induced associated" potential. The latter contributes 
' to the total potential of the disk and thus infiuences the rotation curve. However, 

^ • the nature of the induced associated surface density is geometric rather than ma- 

■ - - ' terial, so it may be perceived as "dark matter" . Applying the derived formulae to 

an exponential disk, we show that we get an asymptotically flat rotation curve and 
a correspondence with MOND, in agreement with observations. 



1 Introduction 

The concept of dark matter (DM), without which mainstream astrophysics and cosmology 
would not be observationally viable, is an acknowledged part of the modern scientific 
worldview. Yet all the observational evidence in favour of its existence is based on 
interpretations of astronomical data coming from distance scales much larger than the 
solar system. In addition, despite its rather flexible nature, the DM proposal faces 
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some real observational challenges. In particular galactic phenomenology, including spiral 
galaxy rotation curve shapes and the seemingly existence of a fixed acceleration scale, 
seems difficult to understand in terms of DM. 

This motivates the alternative approach of trying to explain anomalous galactic ob- 
servations from modified gravity. The most famous of these approaches is the proposal 
known as MOND; i.e., a specific modification of Newtonian dynamics whenever the grav- 
itational acceleration falls below a critical value oq. Fact is that MOND does a very good 
job of modelling rotation curve shapes of spiral galaxies just from their visible matter con- 
tent, using only one freely variable parameter qq. That this feat is possible at all, makes 
the DM approach look suspicious in comparison since no a priori reason is given for why 
a much greater variety of rotation curve shapes is not observed. However, MOND does 
not work so well on larger distance scales such as clusters of galaxies. Besides, relativistic 
extensions of MOND involve arbitrary extra fields reducing their predictive powers and 
thus their advantages vis-a-vis DM. 

Rather than trying to modify gravity by introducing extra fields tailored to fit galactic 
phenomenology, a less contrived approach would be to adopt a general alternative frame- 
work of relativistic gravity and see if galactic phenomena can be correctly predicted from 
first principles. Such an approach can hardly be made to work within the traditional 
framework of metric theories of gravity [1]. However, another possibility is the so-called 
quasi-metric framework (QMF) published some time ago [2, 3]. Quasi-metric gravity is 
rather arcane and it has not yet been shown to be viable. On the other hand, nor is it 
in obvious confiict with observations, even if it would seem so at first glance. The QMF 
is also radical inasmuch as the role of the cosmic expansion is described very differently 
from its counterpart in metric theories of gravity. Specifically, the QMF predicts quan- 
titatively, and from first principles, how the cosmic expansion infiuences gravitationally 
bound systems. In particular, the QMF predicts how the cosmic expansion should affect 
the solar system and that these effects should be observable. A number of observed solar 
system phenomena have been reinterpreted and shown to be in good agreement with 
these predictions as long as no extra theory-dependent assumptions are made [4] . Unfor- 
tunately, these results indicate that general relativity (GR) is fundamentally flawed and 
that interpretations of some solar system observations are crucially theory-dependent. 
Well-known results based on traditional interpretations of such observations may there- 
fore be unreliable. Moreover, questioning traditional interpretations is probably why 
quasi-metric gravity is ignored by the scientific community, despite making a number of 
successful predictions not shared with GR [4]. Besides, if it turns out that quasi-metric 
gravity is able to explain galactic phenomena from first principles as well, this will con- 
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tribute to undercutting the relevance of the DM concept. In that case the legitimacy of 
mainstream astrophysics will be weakened even further. 

The goal of showing the compatibihty of the QMF with galactic phenomena, without 
the assumption of DM, will be partially fulfilled in this paper. That is, we calculate 
the gravitational field of a very thin, rotating disk of dust in the non-relativistic limit of 
quasi-metric gravity. Taking the spinning disk as an idealised model of a spiral galaxy, 
we indicate a solution of some galactic observations without the need of DM. The cor- 
respondence of this solution to that of MOND is also discussed. But first, in the next 
section, a brief introduction to quasi-metric gravity will be presented. 

2 Basic quasi-metric gravity 

The QMF has been described in detail elsewhere [2-4]. Here we include only the bare 
minimum of motivation and general formulae necessary to do the calculations presented 

in later sections. 

The basic motivation for introducing the QMF is the idea that the cosmic expansion 
should be described as a general phenomenon not depending on the causal structure 
associated with any semi-Riemannian manifold. And as we will see in what follows, 
certain properties intrinsic to quasi-metric space-time ensure that this alternative way of 
describing the cosmic expansion is mathematically consistent and fundamentally different 
from its counterpart in GR. 

Briefly the geometrical basis of the QMF consists of a 5-dimensional differentiable 
manifold with topology A^xRi, where Ai = 5xR2 is a Lorentzian space-time manifold, 
Ri and R2 both denote the real line and 5 is a compact 3-dimensional manifold (without 
boundaries). That is, in addition to the usual time dimension and 3 space dimensions, 
there is an extra time dimension represented by the global time function t. The reason for 
introducing this extra time dimension is that by definition, t parameterizes any change in 
the space-time geometry that has to do with the cosmic expansion. By construction, the 
extra time dimension is degenerate to ensure that such changes will have nothing to to 
with causality. Mathematically, to fulfil this property, the manifold A^xRi is equipped 
with two degenerate 5-dimensional (covariant) metrics gt and gf The metric is found 
from field equations as a solution, whereas the "physical" metric gt can be constructed 
from gt in a way described in refs. [2, 3]. (See also section 6.) 

The global time function is unique in the sense that it splits quasi-metric space-time 
into a unique set of 3-dimensional spatial hypersurfaces called fundamental hypersurfaces 
(FHSs). Observers always moving orthogonal to the FHSs are called fundamental ob- 
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servers (FOs). The topology of M. indicates that there also exists a unique "preferred" 
ordinary global time coordinate We use this fact to construct the 4-dimensional 
quasi-metric space-time manifold M by slicing the submanifold determined by the equa- 
tion x° = ct out of the 5-dimensional differentiable manifold. (It is essential that this 
slicing is unique since the two global time coordinates should be physically equivalent; 
the only reason to separate between them is that they are designed to parameterize fun- 
damentally different physical phenomena.) Thus the 5-dimensional degenerate metric 
fields and may be regarded as one-paramctcr families of Lorentzian 4-metrics on H. 
Note that the existence of a "preferred frame" is an intrinsic, global geometric property 
of quasi-metric space-time. (For an isolated system, the preferred frame is the one where 
the system is at rest, see ref. [5].) Furthermore, there exists a set of particular coordinate 
systems especially well adapted to the geometrical structure of quasi-metric space-time, 
the global time coordinate systems (GTCSs). A coordinate system is a GTCS iff the time 
coordinate is related to t via x^ = ct in A/". 

For reasons explained in [2, 3], the form of g^ is restricted. Expressed in a GTCS 
(where the spatial coordinates do not depend on t), the most general form allowed for 
the family gt is represented by the family of line elements (this may be taken as a 
definition) 

, _ _ _ t _ f2 _ 

dst = [N'N, - N^]{dxy + 2-Nidx'dx^ + -^N^Sikdx'dxK (1) 

^0 ^0 

Here to is some arbitrary reference epoch setting the scale of the spatial coordinates, Nt 
is the family of lapse functions of the FOs and ^N'^ are the components of the shift 
vector family of the FOs in (jV, gt). Moreover Sn-dx'^dx^ is the metric of (with radius 
equal to cto) and N—N'^SikN^ . Note that this form of gt is strictly preserved only under 
coordinate transformations between GTCSs where the spatial coordinates do not depend 
on t (but some exceptions to this rule exist). Also note that, since the 5-dimensional 
metrics are degenerate, there are no components of lapse and shift in the ^-direction (i.e., 
there is no motion and proper time does not elapse along the ^-direction) . 

Due to the restricted form (1) of gt, a full coupling to space-time curvature of the 
active stress-energy tensor Tt should not be expected to exist. However, it turns out that 
a subset of the Einstein field equations can be tailored to gt, so that a partial coupling 
exists [2, 3]. The field equations then read (expressed in a GTCS) 

2i?(t)ll = 2(c-2a>|^ + c-V.a>-^W.fc^S + "fn,i?*) 

= ^^(^11 + %), (2) 
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^{t)j± — ^lt)j\i ~ Kt,j — f^T^t)jl- (3) 

Here Rt is the Ricci tensor family corresponding to the metric family gj and the symbol 
'±' denotes a scalar product with — lif, that is the negative unit normal vector field family 
of the FHSs. Moreover, denotes the Lie derivative in the direction normal to the 
FHSs, Kj denotes the extrinsic curvature tensor family (with trace Kt) of the FHSs, 
a "hat" denotes an object projected into the FHSs and the symbol '|' denotes spatial 
covariant derivation. Finally K=87rG/c^, where the value of G is by convention chosen as 
that measured in a (hypothetical) local gravitational experiment in an empty universe at 
epoch to- 

An explicit coordinate expression for may be calculated from equation (1). This 
expression reads (using a GTCS) 

t - - /N ~ t \ - 

where h(t)ij are the components (in a GTCS) of the spatial metric family intrinsic to 
the FHSs. 



3 Axisymmetric, metrically stationary, flat systems 
3.1 Real weak- field solution 



The most general form of gt for an isolated, metrically stationary, axially symmetric 
system can be found from equation (1). Introducing a spherically symmetric GTCS 
{x^, p,9,(l)}, where p is an isotropic radial coordinate, Nt and iV^ do not depend on 
and equation (1) takes the form [5] 



ds, = B 



- (1 - V^p''sm^9){dxy + 2-Vphm''ed<f>dx'' + ^ 



dp" 



+ p'^dn- 



(6) 



where dn'^=de'^ + sin^dd(t)'\ Eo=cto, B=Nf and V=-^^^^. The field equations (2), (3) 
exterior to the source then read [5] 

-"^ 1 ^ 2.. 3p\^ cot^^ 

^0 



P 



p- 



Bp'sm'e\{l-^){V,,f + \{V,or 

"0 P 



(7) 
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3cot^ + % 


S - 



\v,e=Q. (8) 



For flat systems, it is convenient to switch to a cylindrical GTCS 0), where 

^=psin^, z=pcos9. Then equation (6) becomes 



+ 



:i - ^)dC^ + (1 - ^)dz' + 2^dCdz 

^0 



- (1 - V^e){dxy + 2-Ved4>dx'' 



^0 



"0 



and the field equations (7), (8) read 

= Se \{l - )^ + (1 - ^){V,. f - ^4^V,, V,. 



(9) 



(10) 



"0 



+ 



"0 



"0 



1 - 



"0 



B 



El B 



5z 

^0 B ^0 



V,,^0. (11) 



We now assume that the gravitational field is weak, so that we may set V — 0. This 
imphes that equation (11) becomes vacuous and that the right hand side of equation (10) 
vanishes. To find B{^, z), we are thus left to solve equation (10) with the right hand side 
equal to zero. Note that this is equivalent to solving the Laplace equation on a subset of 
the 3-sphere. 

Unfortunately, equation (10) is a non-separable partial differential equation (PDE). 
This means that there is not much hope of finding exact solutions. However, one may try 
to find approximate series solutions for \z\ << Sq. The series expansions in terms of |^| 
should then take the same form as for the corresponding Newtonian problem, recovered 
by letting Sq— J-oo in equation (10). One then gets a separable PDE in the Newtonian 
potential. That problem was solved many years ago [6, 7], yielding a continuous spectrum 
of solutions ^ki^,z) oc Jo{k^)cxp{~'k\z\), where Jo{k^) is a Bessel function of the first 
kind. We are thus led, via correspondence with the Newtonian case, to try solutions 
B{^, z) = 1 + 4$(C) z) built from mode solutions of the form 



MC, z) = '^kiO (l -k\z\ + ^a^ikzf + •••), 



[12) 
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where the ak are constants (but higher order coefficients will in general depend on 
The reason why the ak do not depend on ^, is that any deviation from Euclidean space for 
small 1 2; I occurs at order \z\^ and higher. That is, by integrating the spatial line element 
a distance \z\ « So in the 2;-direction we find (for B^l) 

= z(l + ^—^^ + ■■■). (13) 

V 6(1-|)^§ / 

We now insert equation (12) into equation (10) (with the right hand side set to zero) and 
collect terms to get separate equations for each power oi \z\. To lowest order, the terms 
independent of j^l yield an equation for $fc(^), i.e. 

(1 - ^)^k:^^ + 1(1 - %)^k:^ +akk^^k = 0. (14) 

Similarly, collecting terms of first order in \z\ yields an equation determining the third 
order coefficient of the series expansion in equation (12), and so on for each order in \z\. 

Now, since k is interpreted as a wavenumber on the 3-sphere, it must have a minimum 
value ko = corresponding to the maximum value ^max = Sq. This indicates that, rather 
than the continuous spectrum of solutions found for the Newtonian case, the solution of 
equation (14) should involve a discrete spectrum of solutions $n(0=^fcn(0- ^^^t the 
general solution of equation (14) is (with an=(yk„) 

C-'$n(0 = -CnPnH - C^QM, U = (J^O^S^Tl - l)/2, U=l - 2^, (15) 

where C„, are (dimensionless) constants and Pn{u), Qni^) are Legendre functions of 
the first and second kind, respectively. (Notice that the mode solutions C^(5n('w) diverge 
logarithmically when li— )-±l.) From equation (15) we see that we get a discrete spectrum 
of solutions if n is required to be a non-negative integer. Moreover, we require that k'^i 
should be an integer multiple of /cq, and that in the continuum limit lim„_^ooO;n = 1, so 
we must choose 

1 2n + 1 

Q;n = l-T^> ^ kn^^^ — , n = 0,1,2,..., (16) 

and this choice will also be consistent with the indicated value of k^. Note that, to get the 
correspondence with the Newtonian case, one may choose C]^ = and C„ = e~*^ where 
s is a non-negative real number. Summing over n and using the generating function 

g{u, s) = exp{su)Jo{sVT^) = V -^Pn{u)s\ (17) 

n=0 
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for Legendre polynomials [8] , we get the continuous spectrum of solutions 

'2 



1 L ^ / 

^s{u) (X exp(-s) V -^Pn{u)s^ = exp f - 28^^) Jq hs-^x h - • (18) 

Setting k=^ and then taking the limit S^oo in equation (18), we get back the Newto- 
nian case. 

To find the mode surface densities S„(?i) and T,\{u) corresponding to the specific mode 
solutions —CnPn{u) and — C^Q„(-u), respectively, wc lay a Gauss surface around the disk 
and use Gauss' theorem across it. This procedure is exactly similar to the Newtonian 
case treated in [7]. Assuming a weak field (i?„~l) we find 

^"^''^ ^ 2^^^^ + ^^^"^^^' ^^^^^ = 2^^2^+^^^"^^^- ^^^^ 
Due to the fact that the set of Legendre polynomials Pn{u), nG(— 1, 1), is complete and 
orthogonal [9], it is possible to expand any real surface density S(m) in terms of the 
mode surface densities S„(m). Setting Cn=CSn (where C is some nonzero constant) and 
summing over n, we find 

oo ^ -1 

E(u) = VEn(ii) = ^-^y2-{2n+l)SnPn(u). (20) 

n=0 ^ n=0 

But this means that S(m) is expressed as a Legendre Fourier series [9], and that its inverse 
Sn is the finite Legendre transform [9] of '^^^° I](m), i.e. 

Sn = ^ T.{u')Pn{u')du\ (21) 

SO that 

S(m) = - 5^(2n + l)Pn{u) / E(m')^«("Oc?^^'- (22) 

n=0 "^-1 

The real solution $(w) corresponding to the given surface density E(w) is then obtained 
by summing up the mode solutions — C„P„(ii), i.e. 

oo oo „x 

$(«) - -C^C J2 SnPn{u) = -TtGEo Pn{u) / E(«')Pn(«')c?«'- (23) 

71=0 11=0 

The corresponding series solution B^^^i{u,z) is then found by combining equations (12), 
(16), (22) and (23). The result is 

Br^3\{u,z) = 1 + ^$(h) + ^|^||](k) 



oo 



. z^y^n{n+l)Pr,{u) E{u')Pn{uyu' + 



(24) 
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The circular speed Wreai due to the solution (23) as function of radius of the disk can now 
be found from the usual non-relativistic formula (this is justified for a weak field Si.eai~l, 
see [4]). Expressed as a function of u we find 

n=0 -^-1 

= -r-^J2(n + l)(uPM-Pn+i{u)) E{u')Pr,{u')du'. (25) 

^ + n=0 -^-1 

It cannot be expected that rotation curves found from equation (25) should deviate 
significantly from their Newtonian counterparts (they don't). Thus it would seem that 
the need for galactic DM is the same as for the standard model. However, here we have 
not taken into account possible solutions constructed from the specific mode solutions 
—C\Qn{u). As wc shall sec later, if we do this, new possibilities of getting rid of galactic 
DM open up. But first, in the next section, we need to consider restrictions coming from 
boundary conditions. 

3.2 Boundary conditions 

So far we have implicitly assumed that there are no particular preferences regarding 
the form of E(ii) as long as it is physically reasonable and the resulting ^{u) is small 
everywhere. But is this true? To answer that question, we first notice from equation (23) 
that the real potential at the center of the disk is given by 

$(1) = -vrGHo f: nu')P^W)du' = J'^ (26) 

where the last expression follows from the generating function [8] 

y PMs"^ = , -1 < s < 1, (27) 



m=0 



for the borderline case s = 1. Moreover, equation (26) may be written in the form 



m-c'^, E,.^y_^-^, E,.— , (28) 

where E+ is a weighted average, and where the constant E* sets a specific surface density 
scale depending on the finite size of space. Besides, a second weighted average surface 
density E_, related to the total (active) mass Mt^ of the disk, can be defined by 



1 rEK)£^^ (29) 
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where $(—1) is found from an expression similar to equation (26) by using equation (27) 
for the border hne case s = — 1. We note that E* is a purely geometric quantity, whereas 
E+ and E_ depend on the real surface density profile. 

Next we note that •\/|$(— 1)| represents a specific (non- vanishing) velocity scale. This 
indicates the possibility of defining some quantity ;^[a-(— 1) — a-(l)]~E* with the property 
that it relates $(1) to ^|$(-1)| via a definition similar to equation (28). By combining 
equations (28) and (29) such a relationship may readily be found. That is, we may define 




*(l) = -c,/|*(-l)&^.-c^¥FT^^Mzl^5Wl. (30) 



Notice that equation (30) involves the geometric quantity E*, but since $(1) = ^$(— 1), 
the analogous relationship between $(1) and $(—1) does not. Also notice that the factor 
7r/-\/2 is included into the definition (30) since a'(— 1) — a(l) should be more similar to a 
mode surface density (see equation (19)) rather than to a weighted average like E_|_. 

A definition similar to equation (30) may be made for the contribution $>u(l) to $(1) 
from the part of the disk interior to some arbitrary coordinate u. The purpose of such a 
definition is to construct a new "associated" surface density a{u). That is, we may define 

-c/W^^Mh)^, (31) 
or equivalently (where the constant a'(l) must be determined separately, see below) 



/X /X 1 /S, f^T.iu')du' , , , , \/2 /E, - 

We see from the definition (32) that a{u) is increasing from the center of the disk and 
outwards. Thus, a{u) could be interpreted as some kind of "inverted" surface density. 
This means that (y{u) should not be considered as an independent, gravitating source. 
Rather, (y{u) should give some restrictions on the possible forms of S(m). Such restrictions 
can be found by requiring that (j{u) should be linearly related to E(m), so that the 
corresponding potential can be written as a linear combination of $(m) and some constant 
potential. Then a{u) is not independent. A "preferred" form of ^{u) can thus be found 
by requiring that ^{u) — + AE(m) — a"(— 1), where A is a constant, or equivalently 
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a(M) = a(-l) - ^— ^ S(m) -S(-l) , a(l) = E(-l)^— ^. (33) 



E(l) • 



(We see that in this case, since E(— 1) should be neghgible, a{l) should be also.) Equation 
(33) is then an integral equation determining the most basic form of the real surface 
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density for an isolated disk. To find exactly what this form is, it is convenient to 
turn equation (33) into a first order separable differential equation by taking the derivative 
w.r.t. u at both sides of it. Solving this equation is straightforward, and the result is an 
exponential disk, i.e.. 



E(m) — E(l)exp 



(S^ - ^) V(l - ^)/2] =E(l)exp [ - I V(l - «)/2] , (34) 



where ,^d=2oS+/[S(l) — S(— 1)] is the disk length (at epoch to)- This result answers 
the question we posed at the beginning of this section; the simple requirement that S(m) 
and a{u) should be linearly related implies that there is a particular preference regarding 
the form of S(m). That is, it would seem that the exponential disk should represent 
a preferred surface density profile among all the possibilities that might exist. This is 
confirmed observationally, since an exponential surface density is the hallmark density 
profile of the outer regions of spiral galaxies. We will return to the exponential disk in 
section 4. 



3.3 The induced solution 

Contrary to the Legendre polynomials Pn{u), the functions Qn{u) are not polynomials, 
and nor do they constitute an orthogonal set for «€(— 1, 1). Rather, the functions Qn{u) 
can be separated into two subsets depending on whether n is even or odd. That is, each 
function with odd n is orthogonal to every function with even n and vice versa. On the 
other hand, functions within each subset are linearly dependent. This can be easily seen 
from the formulae [81 



Qn{u)Qn{u)du 



6 ^ Z^fe=l fe2 



1 / 2n + l 

Qn{u)Qm{u)du = — - — — tt( y'-y] It, T^T^m. (35) 

i (m — n)(m + n + 1) \ ^ 

The problem now is to construct a solution $'(«) from the mode solutions —C^Qniu) such 
that $(w) and ^^{u) are linearly independent, i.e., we require that J_^<^'{u)^{u)du^O. 
However, since there are only two sets of linearly independent mode solutions, it must be 
possible to find many such solutions by summing over different numbers of mode solutions 
(in general, at least two mode solutions must be included, one from each linearly inde- 
pendent set). This means that, merely requiring linear independence is not sufficient to 
arrive at a unique solution ^\u). However, a unique linearly independent solution $'(u) 
can indeed be found by summing over all the mode solutions. We will call this solution 
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the induced solution, since it is found indirectly by summing up all the mode solutions 
C!^Qn{u) such that every term in the mode sum has a linearly independent counterpart 
CnPn{u) from equation (23). Moreover, the requirement of linear independence is not 
trivial since it forces the constants C\ to be dependent on Sn and thus the real surface 
density. The corresponding surface density S^(ii) will be called the induced density. The 
induced density is not real, but could still have physical consequences indirectly. 

The induced solution $'(«) obtained by summing over all the mode solutions can be 
found by assuming that the coefficients C,\ can be written in the form C\ — C^Sn, where 
C" is a normalisation constant. Then, using the formula [8] 

1 2 ( j^'jn+m 

Pn{u)Qm{u)du = ^ ■ — , Tly^m, (36) 

_i [n — m)[m + n + 1) 

it is easy to see that 

fl oo oo (_l)n+m 

J-1 fr^^o (n-m)(m + n + l) 

from the obvious antisymmetry obtained by permuting the summation indices. To com- 
pletely specify the solution $^(m), it remains to specify the normalisation constant C\ 
But the only natural choice is really to set — C. We thus have that 

Q = C„ = CSr, = f i:{u')Pn{u')du'. (38) 

The solution ^^{u) is now completely specified, and we have 

oo QO ^\ 

^\u) - -c'C J2 SnQn{u) = -ttGSo / n^')Pn{u')du' . (39) 

n=0 n=0"^~l 

Equation (39) may be written in a more convenient form by using the identity [8] (the 
integral being defined by its Cauchy principal value) 

(40) 

so that by using equation (23) we find that (again using Cauchy principal values) 

Prom equations (19), (22), (38) and (40) we can also find the induced density S'(xi) 
corresponding to the induced solution ^^{u), i.e., 

^'N = E ^(2- + mn{u) = \ f (42) 
^ 27rG=,o 2 J_i u-s 
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However, there is a fundamental problem with the induced quantities. That is, both 
and S'(ti) will in general contain generic logarithmic divergences (typically located 
at u— ±1) inherited from the mode solutions Qn{u). This means that ^^{u) and T,^{u) 
cannot be used directly as physical quantities. But S'(ii) may be used indirectly, since 
it is possible to construct a new, geometric, singularity-free surface density from it. We 
will show how to do this in the next section. 



3.4 The induced associated potential 

Although S'(m) cannot be used directly, it still may have physical significance. The reason 
for this is that it is possible to construct a new surface density d-\u) from S'(m) using the 
results found in section 3.2. There, the definition (32) of a{u) was motivated from the 
possibility of restricting possible forms of S(m) due to boundary conditions. However, 
since a{u) is by construction unphysical, there never was any reason to interpret it as an 
independent gravitating source. On the other hand, there is always the possibility that a 
definition similar to (32), but with S'(ii') substituted for S(u'), may have the properties 
that makes it possible to use it as a gravitating, geometric source. That is, we may 
construct the so-called induced associated surface density a'^(w), defined by 



Note that, by integrating over S'(m'), its generic logarithmic divergences disappear, so 
that d'\u) does not contain such singularities. Moreover, since it is possible that £'(ii) 
may change sign somewhere in the interval u&{—l,l), d'^u) may take a form similar to 
some physical surface density profile. Finally, unlike a{u), a'^{u) can certainly not be 
algebraically related to E(ii). Thus d'^{u) may be considered as a physical, independent 
gravitating geometric quantity, but not as a material density. 

We may now use d'\u) as a source in equation (23) to get a new potential, i.e., the 
so-called induced associated potential ^(li) defined by 



n=0 -^-1 



(44) 



Since d-^u) does not contain generic divergences, we see from equation (44) that '^{u) 
should be non-singular everywhere, so it may be accepted as a physical quantity. Thus 
to any real surface density profile S(u) and its corresponding real potential ^{u), there is 
always associated a surface density d-^u) (playing the role of galactic "dark matter"), and 
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its corresponding potential '^/{u). These quantities should always be considered together 
when making predictions. 

The series solution B{u, z) containing both real and induced associated contributions 
is then given by 

- ^z' V n{n + 1)PM f inn') + o\u')) P„iu')du' + ■■-. (45) 

Note that the second order term in equation (45) is divergent for w = 1, Zy^O. See the next 
section for more comments. Moreover, the (non-relativistic) circular speed Wdrc calculated 
from equation (45) is due to both real matter and induced associated "phantom" matter, 
and similar to equation (25) we find that 

wlM = f^n + 1) (uPniu) - Pn+liu)) j'^ (s(«') + a' (u)) P^iu' ) du' . (46) 

We see from equations (41) and (42) that the induced solutions and densities are obtained 
by integrating real solutions and densities over the whole disk; thus these quantities 
describe non-local, collective properties of the system. This is also true for a\u). That is 
why any extra gravitational acceleration obtained from ^(u) should not be interpreted 
as a fundamental modification of the Newtonian force law such as in MOND - rather the 
extra acceleration should be seen as an emergent property of the whole system. 



3.5 An important transformation 

Given as input the real surface density S('u), one should now in principle be able to 
calculate the real potential $(«), the real scries solution Bj.cs.i{u, z), the corresponding 
induced quantities S'('u) and $'(m), the induced associated quantities (t'('u) and ^(m), the 
total series solution B{u, z), and finally the rotation curve from equation (46). However, 
as seen from equations (23) and (44), the expressions for $(m) and ^(m) contain an 
infinite sum over (a product of) Legendre polynomials. Since this sum will in general 
converge slowly, its presence makes numerical calculations quite awkward. Fortunately, 
it is possible to rewrite this infinite sum in terms of an elliptic integral, making numerical 
calculations much easier. The key to this important transformation is using the generating 
function [10] (|s| < 1) 

J2 Pm{u)Pm{u')s^ , "^"^ (47) 

"»=0 ^ ° W 1 - 2s (uu' + yJ{l-u''^){l-U^)cOSU?j + S2 
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for the borderline case s = 1. (Note that, for the special case u' — 1, we get back equation 
(27).) 

We now use equation (47) to rewrite equations (23) and (44). Interchanging the sum 
and the integral in these equations, adding them and then using equation (47), we get 



$(m) + -^{u) 



Gl 



E(m') + a\u') 



du'duj 



^ -^0 J-^ - uu' - 7(1 - ii'2)(l -^cosa; 
E(m') + aHu' 



-V2GE0 



K 



2^(l-«2)(l-^/2) 



1— !iM'+7(l-«2)(l 



-1 



t{u') + a\u')\K(^^^^?^^)du' 

\ ' ' V / \ i—uu'+\u—u'\ J 



^1 7(1 -m'2)(1 -m2) + ^ I - uu' - s/ {I - U''^){1 - U^) 



m 



where K{k)= J^^^ dO/y/l — k'^sin'^O is the complete elliptic integral of the first kind [8]. 
The last form (see [10]) of equation (48) is a httle better to use for numerical purposes. 

Similarly, it may also be tempting to interchange the infinite sum and the integral 
present in the second-order term of equation (45), and then transform the infinite sum 
into an integral. However, as we shall see, this procedure does not work for higher-order 
terms, since all will be divergent. To illustrate this, a somewhat lengthy calculation using 
equation (47) yields 



J2 m{m + l)P^{u)Pm{u') = lim J] Pm{u)P^{u')s 



m=0 



E{k) - (1 - uu' - ^{l-u^){l-u'^))K{k) 



k^=- 



d ,d 



ds ds 



m=0 



2v/(l-«2)(l- 



w 



I2\ 



2^f2^l\u - u'\^l - uu' - ^1 - u'''){l - u^) I - uu' + ^ [l - u''){l - u'^) 



:,(49) 



where E{k)= ^'^ ^/l — k'^s\v?6d6 is the complete elliptic integral of the second kind [8] . 
To get an expression for B{u,z) (for small |2;|) more suitable for numerical calculations, 
one may now insert equation (48) into equation (45). However, if one tries to insert 
equation (49) into the quadratic term, interchanging the infinite sum and the integral, 
it is straightforward to see that this term will diverge. This behaviour is quite similar 
to the Newtonian case (where the potential takes the form of a double improper integral 
[7]) if one tries to expand the Newtonian potential in powers of interchanging the two 
improper integrals implies that all terms with power >2 will diverge, yet if all the terms 
are included, the resulting exact result is finite. But only including the linear term in 
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the series expansion approximates the exact result well for small enough \z\. One expects 
that this is valid for equation (45) also, so that skipping terms with power > 2 is justified 
for small l^j. Note that in equation (45) (also just as for its Newtonian counterpart), all 
terms of order >2 will diverge on the z-axis for z^Q. 

Of course, if \z\ is not small enough, the second order term (and possibly higher-order 
terms) must be calculated in some admissible way and included. That is, it would be 
preferable to find some other way of rewriting the series expansion in \z\ (i.e., without 
using equation (49)), such that all terms can be expressed in a form suitable for numerical 
calculations. Fortunately, such a form for the second order term can be found from 
equations (45) and (10) (with V — Q). From these equations we find that (for small |2;|) 



B izz — 



"0 S "0 

;i + «)($,„„ +^,uu ) - 2m($,„ +^,u . 



0{\z\l 



(50) 



so that the first few terms of the rewritten series expansion read 



B{u,z) = l + -($(ii) + *(xi)) 

4^2 



+ 



At:G, 



\{i:{u)-ra'{u)) 



(1 + «)($,„„ +'^,uu ) - 2m($,„ ) 



+ 0(\z^ 



(51) 



Note that the second order and higher-order terms are still expected to diverge in the 
limit rt— )-l, 2:7^0. This means that, to do calculations to a given accuracy at some fixed 
2:7^0, one needs to include ever more higher-order terms into the series expansion when 
moving towards the z-axis. 

Finally, interchanging the sum and the integral in equation (46), yields, after some 
tedious calculations, that 



u 



l + u 



T.{u') + a'{u') 



K 



2i/(l-u2)(l_„'2) 



1-mm' + Y^(1-«2)(1-„'2) 



du' 



l + u 



I [u — u 



t.{u') + aHu') 



E 



2^(l-«2)(l-«'2) 



l-uu'+y(l-u2)(l_u'2) 



du' 



-u'\^l-uu' - ^J{l-u''^){l-u^) 



(52) 



\u 



Inserting any given surface density S(-u) and its induced associated surface density cr'(ti) 
into equation (52) now yields the full rotation curve of the disk. 
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4 The exponential disk 



4.1 Approximate solutions 

For spiral galaxies, the observed general trend is that surface brightness (and thus the 
luminosity due to stars in the disk) falls off exponentially from the center and outwards. 
Therefore, to explain spiral galaxy rotational curves without dark matter, one is required 
to assume that surface density profiles of stars are proportional to luminosity profiles. 
This means that, in the general, but idealised case of a truncated disk of extension ^q, 
we must have that 



Ecexp( 




^ f Eeexp( - f l>u>Uo, 







U < Uo, 



(53) 



where Ec is the central surface density and is the disk length (both taken at epoch 
to)- Note that, since ^d^^o, S(^) falls off so fast that we are justified in treating spiral 
galaxies as truncated disks for sufficiently large ^o, i-e., for ^o^Cd- For computational 
purposes, serious errors are hardly made if we choose to lie somewhere in the interval 
20 Cd^Co^^^o. On the other hand, truncation inevitably yields truncation singularities. 
However, their contributions to calculations are usually negligible. In what follows, we 
will assume that the disk is not truncated, but the results are not significantly affected 
if it is. 

Now we may calculate the induced surface density E'(^) associated with the exponen- 
tial disk. To do that, we compute the integral 



(u-s) ^exp( - =^ 

^ ^d 



/So 1-u 



-p(-|/¥){-[-|( 



1-u 



)] 



Ei 



Ei 



— 

-^d 



1 - s 



1-u 



ds 



]} 



1-u 



— exp 



/So 1 



Ei 



So /I 
^d 



u 



— exp 



Ei 



u 



where the approximation holds as long as u is not very close to —1, and if ^d^^Q. 
Ei(a;) is the exponential integral defined by [8] 



(54) 



Here 



EiCa;)^ - 



exp(— s) 



ds 



f 

J — oc 



exp(s) 



ds, Ei(— a;)= — 



exp(— s) 



ds, X > 0. (55) 
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We then find from equation (42) that 

S^(0«^ - Y{exp(e/ed)Ei(-e/ed) + exp(-e/ed)Ei(e/ed)}, < ^ < Sq. (56) 

Note that for (moderately) large distances, unlike ^{^), ^^{0 tends slowly towards zero 
(from below) with increasing ^. Besides, note that S'(,^) is positive for small ^ and that 
it diverges logarithmically towards the origin. (As can be seen from equation (54), the 
exact expression for E'(0 also has a similar singularity when ^^Sq.) 

The induced associated surface density a-^(^) is now straightforwardly found from 
equations (43) and (56). Neglecting the very small quantity o-'(O), we find that (setting 
Mt(jf=:i27rEcCd' i-®-' neglecting corrections coming from the finite size of space) 



^{exp(-e/ed)Ei(e/ed) - exp(e/ed)Ei(-e/ed)}, 0<^<^o. (57) 



Note that d-^{^) (even using the exact expression for S^(0 found from equation (54)) is 
non-singular and non-negative everywhere. Thus ct'(^) may formally play the role of DM. 

Next we find an approximative expression (valid if ^<CSo) for the real potential $(0 
from equation (48) (omitting a"'(u)). Splitting the integral up into two improper parts, 
we find that 



$(0~ - 4GEeed 



a;exp( 



-x)K 2 



dx 



X + e/Cd 



(58) 



Note that, as can be readily verified numerically, $(^) as calculated from equation (58) is 
essentially identical to its Newtonian counterpart $n(0 (valid for an infinite exponential 
disk), given by [7] 



(59) 



where I,^{x) and K,^{x) are modified Bessel functions [8]. Similarly, the corresponding 
circular speed Wreai can be found approximately from equation (52) (omitting d'\u)), and 
this is also numerically very close to its Newtonian counterpart. That is, again splitting 
up integrals into improper parts, we find that for ^-CSq, 



^Li(0«^2GEeed 



31 

ed 



xexp{—x) 



K 2 



e(2^§§^ 



X + e/Cd 



X - e/Cd 



dx. 



(60) 



Note that the last term in equation (60) diverges for each integral, but that added to- 
gether, the sum converges (taking its Cauchy principal value). 
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An approximate expression for the induced associated potential is given by 



TT" V =^0 



for ^<^So. Note that the magnitude of ^{^) depends critically on the upper hmit of 
integration (if ^~So, a more accurate expression for ^(^) may be found from equation 
(48)). Also note that ^(^) is non-singular everywhere, and that if Ec is small enough, 
1^(^)1 may in principle dominate over |$(^)| for the whole disk. Finally, the full rotational 
curve may be found approximately from the formula (^<^So) 



"^circV 



iO^y^reM) + + ]^{exp(-x)Ei(x) - exp(x)Ei(-x)} 

Id 



Numerically, it is found that with increasing ^, the integrals converge rather quickly 
towards a constant factor of 27r, so for ^ above about 8-10 disk lengths the induced 
associated contribution changes very little. This means that, for large distances, the 
expression for Wcirc(0 does not depend significantly on the upper limit of integration, 
and that one automatically gets an asymptotically flat rotation curve. Also notice that 
for small distances, the integrals yield a negative contribution to WcirdO- 

4.2 The correspondence with MOND 

The most basic feature of MOND is the postulation of a fundamental acceleration scale oq, 
observationally estimated to ao~(l — 2)xl0~^°m/s^. That is, for Keplerian accelerations 
well below oq, the Newtonian acceleration on = in a spherically symmetric system 
is replaced by the MOND acceleration om = VGA^ao _ the help of results derived in 

the preceding section we may now easily find a correspondence with MOND and actually 
calculate qq in the case of an exponential disk. Prom equation (62) we find that the 
asymptotic rotational speed and the corresponding Keplerian acceleration are given by 




^ a{0^-\— — T, C>Cd- (63) 

We are now able to compare equation (63) to the MOND acceleration. We then get a 
correspondence if 

2c^ 2 

ao = -—- = — c//o~l■5xl0-^°m/s^ (64) 
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where the expression is vahd for an exponential disk at epoch to (with Hq as the cor- 
responding Hubble parameter), and where the final, numerical result is valid for the 
present epoch. (Note that qq is predicted to decrease with time, since the Hubble param- 
eter decreases as the inverse of cosmic epoch.) This means that there is a correspondence 
between the model of a thin disk presented in this paper and MOND for the asymptoti- 
cally flat part of the rotation curve. In particular, the fact that MOND successfully flts 
the observed baryonic TuUy-Fisher relation [11] means that quasi-metric gravity does as 
well. That is, from equation (63) we get that 

and since the constant of proportionality in the relation M^yOCW^ is measured to be about 
50 MQS^/km'^ [11], it is straightforward to check that this result is in very good agreement 
with equation (65) (for the present epoch). This justifies the definition of o'(— 1) made in 
equation (32). Note that, what enters into equation (65) is the active gravitational mass 
Mt (at epoch to); this increases linearly with cosmic epoch. Moreover, since the Hubble 
parameter decreases as the inverse of cosmic epoch, this means that Woo is constant; 
i.e., that it does not evolve with epoch. In other words, the quasi- metric model predicts 
that the sizes of metrically stationary galactic disks increase with epoch, but such that 
rotational speeds arc unaffected. 

Moreover, a FO at an arbitrary epoch ti, using his locally measured values and 
l^-f^o of the gravitational "constant" and the Hubble parameter, respectively, is predicted 
to measure the same slope of the local baryonic TuUy-Fisher relation as will a FO at 
epoch to- This means that the slope of the baryonic TuUy-Fisher relation is predicted 
not to depend on epoch, and this seems to be confirmed by observations [12, 13]. (The 
analyses presented in these papers are based on the standard cosmological framework, so 
some inconsistency with the QMF might be expected. But as long as theory- dependent 
gravitational physics (i.e., inconsistent with the QMF) is not used to infer masses, this 
should not matter too much.) 



5 The general weak-field, axisymmetric case 

Having analysed a fiat disk, the question now is if the results of the previous sections 
can be extended to the general, weak-field axisymmetric case. To answer that question, 
we will follow the standard method of breaking up an arbitrary axisymmetric matter 
distribution into a series of concentric, spherical shells (of negligible thickness) with cor- 
responding surface densities. The total potential at a given point will then be the sum of 
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the potentials due to the collection of shells. Except for the limitation to axial symmetry, 
this procedure is the counterpart to the derivation of a general multipole expansion for 
the Newtonian case as given in [7] , which we will follow closely. 

The mathematical task is then to solve equation (7) (with V — 0) interior respectively 
exterior to a given isolated shell, subject to suitable boundary conditions. Since equation 
(7) becomes separable we may write solutions in the form B{p,9) — 1 + 2c~^$(p, 9), were 
c~^$(p, 6)= — F{p)G{9). General mode solutions Fjj±{p), G/j±{9) are given by equations 
(A. 3), (A. 4) in appendix A. We will select specific mode solutions Fn{p), Gn{9), where 
n = 0, 1, 2, ... are whole non-negative numbers, such that 



Pn{cos9), (66) 



where An and Bn are constants. We now construct solutions $int(p, ^), ^ext(p, ^) respec- 
tively interior and exterior to a shell located at p = by summing up suitable mode 
solutions. That is, by requiring that the interior solution should be regular at the center 
of the shell we find that 



C-^$int(p, ^) = - J— E ^-^r"' (V ^ " |2)^n(cOS^), P < Ps- (67) 

V P n=0 2 V "0 

The exterior solution is found by requiring that $ext(So, 9) must vanish. So, for p > ps, 

P„(cose),(68) 



C-^$ext(p,^) = -W-$^ 




n=0 



1 - 1|) + b^^pT^ 



1 - 



'='2 
"0 



where the constants Dn are given by the expression D„ = —BnPi ^ (0)/Pi ^ (0), which 

2 2 

vanishes for all even n (see equation (A.5)). 

Now the surface density Es{ps,9) of the shell can be expressed as a sum of modes; 
this yields counterparts to equations (20)- (22) valid for a thin disk. We thus have 

Ss(Ps, ^) = - X](2n + l)Esn(Ps)Pn(cose), Esn(Ps) = / ^s{Ps, 9') Pn{cos9')sm9' d9' .{69) 

n=0 "^O 

To determine the constants A„ and S„, we now apply Gauss' theorem across the shell. 
Assuming a weak gravitational field (B^l) we then have 



'1- 



"0 



V op / p=ps \ op J 



dp 



P=Ps 



47rGEs(ps, I 



(70) 
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Since the potential must be continuous over the shell, we have that $ext(Ps, d) — ^int(Ps, d). 
Furthermore, due to the orthogonality of the Legendre polynomials, this means that 

2 

A second equation relating the constants An and Bn can be found from equation (70) by 
inserting equation (69). Calculating the derivatives and using recurrence formulae valid 
for Legendre functions, one finds that (since the Legendre polynomials are orthogonal) 




One may find explicit expressions for An and Bn from equations (71) and (72). We get 

An = {2n+l)M^^9n{Ps), Bn = {2n + l)M^^fn{p.), (73) 



/n(Ps)^ 



Pi 



n- 



■^2 ' 
"0 



(74) 



gn{ps 



2 



pr-\y.) pr'Ho) 



fn{ps)- 



(75) 



Equations (73)- (75) may now be inserted into equations (67), (68) to get complete ex- 
pressions for $int(p, ^) and $ext(p, ^)- This yields the potential generated by a single, 
infinitesimally thin shell. 

To evaluate the potential generated by an entire collection of shells filling space, we 
let SJ^siPs, d) and SJ^sniPs) denote the relevant quantities for a shell lying between ps and 
Ps + 5ps. Prom equation (69) we then have (inserting 5T,s(Ps-, — QmiPs, Q')^Psl \ h- ~ ^) 



^'^snips) 



Jo 



gm{Ps,9')Pn{cos9')sm9'd9' 



Sps 



2g„,n{Ps)Sps 



1 PS 
-I- ~ 



1 Ps 

-I- ~ 



(76) 



22 



where QmiPs, d) is the coordinate volume density of (active) mass [4, 5]. Substituting equa- 
tion (76) into the complete expressions for the corresponding potentials 5^int{p,d) and 
(^^ext(p, 0) and integrating over ps, we find that (with Xs=-^, y=-Jl — ^ and Q^^,= ^^Q„i ) 

P So 

Ps=0 ps=p 

^(2n + l) mr/D^+l/ N ^5 D-"-i/ \ 

Z^^- -Pn{cos9) \P, \y) 2-^ P, '{y) 




X / ; . r -t 1 



2 

=0 



(77) 



Equation (77) is the axisymmctric counterpart to the similar general multipole expansion 
formula in Newtonian theory; this is given explicitly in [7]. These two expressions have 
a correspondence in the limit Sq— 7>oo. For example, the potential at p = is given by 
c~^$(0) = —y— Xsgjj^o{xs)dXs, the same as the Newtonian expression in said limit. 
This means that as expected, equation (77) can not describe DM-cffccts. 

Moreover, due to the boundary condition at p = Sq, there is no obvious way to define 
an associated (volume) density as a counterpart to equation (32). But if no associated 
density can be defined, one cannot motivate a definition similar to equation (43) for 
the induced associated (volume) density. So it would seem that, the explanation of DM- 
effects presented in this paper should exclusively be connected to disks or other structures 
that do not satisfy the boundary condition $(So) = 0, answering the question asked at 
the beginning of this section. 

However, DM-effects are seen in galaxies other than spiral galaxies. In particular, ob- 
servations indicate that dwarf spheroidal galaxies are the most DM-dominated systems 
ever found [14] . Yet the general observational status for the existence of DM in ellipti- 
cal galaxies is more complicated than for spiral galaxies, since some ordinary elliptical 
galaxies apparently lack significant amounts of it [15], while others seem to have plenty 
[16]. It still remains the challenging task of explaining these observations without DM. 



6 Gravitational lensing 

For a sufficiently weak gravitational field, rotation curves can be calculated accurately 
enough using the auxiliary metric family gt rather than the full "physical" metric family 
gf However, when calculating deflection of light, or gravitational lensing, it is not suf- 
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ficient to know gt, even if the gravitational field is weak. But to calculate gravitational 
lensing, it is fortunately not necessary to know in full; as we shall see, a suitable 
approximation will be sufficient. 

The general formulae describing the transformation gt— >gt are given by [2,3] 



9it)Q0 = ~ ^) 3{t)oo 



9m 
9m 



9m + tl{i- 



4- 



(78) 
(79) 

(80) 



where 66=^6^^ and e^=j^e\dx^ are unit vector and covector fields, respectively, along 
the 3- vector field bjr found from the set of linear algebraic equations 



h% - 2a^ = 0, 



(81) 



and where v=c~^y^ aj^kf^^bj^iby. For a weak gravitational field and for distances much 
smaller than Sq (at epoch to), we have that f^c, so to a good approximation we may 
neglect terms of order 2 or higher in the small quantity v/c. This means that we may set 
(assuming 7Vj~0) 



9{t)oo^9{t)oo, 



9it)oj^0, 9m^9m + ^72 T^i^r 

Iq c 



(82) 



To get exphcit formulae, it is easiest to solve the set of equations (81) using standard 
methods, in spherical coordinates. One may then transform to cylindrical coordinates. If 
one additionally assumes a weak field in vacuum, and neglects all terms proportional to 
Sq^, one may give bjr and bjr as series expansions in z. Terminating the series after the 



linear term, yields (using the fact that in vacuum, B,^ 
(10)) 



B,^^ —^B,^ from equation 



{B,ii^B,^)B,^^+(B,^,y 
-8,2 —-8,^2 B,^ 



+ -b^^ + 0{z% 



(83) 



(84) 



where the derivatives may be calculated approximately to first order in \z\ from the 
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expression 

2/_.. 47rG' 



B{^,z) « l + -($(0 + *(0) + ^(s(e) + aXe))| 



1 



*,€€+^*,e+*,$€+^^,«]^' + 0(|2;|3), (85) 

obtained from equation (51) and the approximation for B,zz shown above. Furthermore, 
since the quantities e| and may be found to the relevant accuracy from the definition 



- - V' h ■ I 

e? = V)y4=^(to)»i^=^> |b^|=y6^fe6^, (86) 

we may also calculate the relevant approximation of the desired quantities g(t)i3 from 
equation (82), i.e.. 



gm^9m + T^iS-^W^i + ^71\ u us ^nhTr (87) 



From the geodesic equation, we may now calculate the gravitational bending of a light 
ray grazing the plane of the disk, since in this case, l-^l is small enough so that the approx- 
imation given in equation (85) is valid. To deal with the opposite situation, where the 
light path is nearly orthogonal to the disk plane, the approximations given in equations 
(83)- (85) may not be sufficient; then one must include terms of higher order in \z\. 

Finally, we note that in GR, the weak field form of the metric outside the disk is 
assumed to take the form 

ds'' = -(1 + ^$N)(cix°)2 + (1 - |$N)(cix2 + V + dz"), (88) 

where $n is the Newtonian potential for the sum of visible and DM. But in general, the 
quantities g{tQ)ij as found from equation (87) will not correspond to their counterparts 
given by equation (88). This means that any observationally based mapping of DM 
distributions using gravitational lensing, assuming the weak field approximation obtained 
from GR, is explicitly theory-dependent and may give misleading results. 



7 Discussion 

For many years, it has been known that galactic dynamics is incompatible with a straight- 
forward application of Newtonian theory to visible matter. However, the most glaring 
discrepancies between observations and theory can be removed by assuming the exis- 
tence of galactic DM. Since the introduction of DM can be done without making radical 
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changes to the standard theoretical framework underlying mainstream astrophysics, this 
is currently the preferred approach. On the other hand, MOND interpreted as an em- 
pirical recipe, has an impressive successful record when predicting galactic phenomena. 
But the connection between MOND and fundamental physics has been unclear so far. 

Contrary to other approaches, the explanation of some galactic phenomena given in 
this paper has not assumed any empirical aspect of galactic dynamics as input to the 
model. Rather, while formally belonging to the "modified gravity" category, the model 
comes directly from the weak field approximation of the QMF, without any extra modi- 
fications of the theory. (The only extra assumption made, is that the induced associated 
surface density a^u) should be treated as a gravitating source in the field equations.) 
The main reason why this is possible is that according to the QMF, the Universe is finite 
and "small", so that boundary conditions depend crucially on the shape of the matter 
distribution. (A finite and "small" Universe is incompatible with cosmological data as 
interpreted within the standard framework, therefore the DM explanation given here 
is compatible with the weak field limit of the QMF but not with the weak field limit 
of standard cosmology.) In particular, for a fiat disk we found that = —1)7^0 in 
equation (29), defining a specific velocity scale dependent on the total mass of the disk. 
To be able to define the associated surface density o'(m) in equation (32), it is essential 
that this velocity scale docs not vanish. However, by construction, it docs vanish for the 
general axisymmetric matter distribution considered in section 5, so fiat disks seem to be 
an exceptional case. Anyway, that case should apply to all thin disks (even if they are 
not exactly flat). 

The other crucial feature is the existence of an induced matter surface density S'(m) 
directly dependent on the real matter surface density as shown in equation (42). This, 
together with the existence of a{u), is sufficient to define the induced associated surface 
density d-^u) playing the role of DM. The introduction of d-^u) may seem "contrived" to 
some people. However, the facts that the induced associated surface density correspond- 
ing to an exponential disk automatically yields an asymptotically flat rotation curve and 
a correspondence with MOND are calculated results, and not put in by hand a priori. It 
was not at all obvious that these results would be possible. 

It has been claimed [17] that observations of gravitational lensing in the colliding 
clusters 1E0657-56 (the Bullet cluster) represent a "direct" detection of DM, since these 
observations indicate that the DM is associated with the regions containing the held 
galaxies rather than with the regions containing the more massive gas making up the 
bulk of the cluster. However, as we have seen, in the QMF, the existence of any sort of 
"phantom" matter density similar to d'\u), playing the role of DM, is crucially dependent 
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on the shape of the matter distribution. This means that, e.g., a large nearly spherical 
or spheroidal mass distribution of gas, should not necessarily be associated with much 
DM as inferred from gravitational lensing. So, since the colliding gas clouds in the 
Bullet cluster shown in [17] do not seem to have shapes that could in any way resemble 
disks, this might be a natural explanation of why they do not seem to be associated 
with much DM. But of course, further justification of this explanation will be necessary, 
together with an explanation of why dwarf spheroidals and some elliptic galaxies seem 
to be DM-dominatcd. Anyway, the mere existence of such an explanation shows that 
the interpretation of the Bullet cluster observations is not theory-independent, so citing 
them as definite evidence of the existence of DM is unjustified. 

In light of the results found in this paper, it seems that some lines of argument 
favouring DM over modified gravity have been shown to be invalid. First, there now 
exists a natural correspondence between MOND and fundamental physics. This indicates 
that at least some galactic phenomenology has its basis in geometry rather than in the 
properties of some unknown exotic particle. Second, while the correspondence with 
MOND works for spiral galaxies, this does not imply that such a correspondence is 
necessarily valid for other types of galaxies or galaxy clusters. This means that it may 
be possible to share MOND's successes but not necessarily its failures. (Further work 
should be done to see if this is indeed the case.) Third, some specific observations (of, 
e.g., the Bullet cluster) have been hailed as the ultimate "direct proof" that DM prevails 
over modified gravity. But since standard interpretations of such observations are model- 
dependent, this is not true. 

There is now near consensus that galactic phenomenology has it basis in some suitable 
nonbaryonic particle, despite the fact that no such particle has ever been detected directly. 
Furthermore, many astrophysicists are of the opinion that this is not necessary, i.e., that 
indirect observations are sufficient to rule out any alternative explanations. This means 
that for many astrophysicists, taking the focus off DM and recognising the merits of 
modified gravity, is a radical move that is out of the question, regardless of how contrived 
and ad hoc explanations based on DM might be. This is a perilous attitude, since 
interpretations of indirect observations are often crucially theory-dependent. Another 
example of this is interpretations of some indirect observations made in the solar system 
that are also presented as indisputable facts, even if alternative interpretations are not 
ruled out [4]. So it would seem that, for at least some parts of astrophysics, there is a 
serious lack of critical assessment of basic assumptions underlying mainstream knowledge. 
One may hope that this will change in the future. 
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A Boundary behaviour of mode solutions 

In this Appendix, we list mode solutions of the vacuum field equations for axially sym- 
metric, metrically static, isolated sources with various "pure" values of their multipole 
moments. As we shall see, these solutions can be classified into two groups; those that 
admit the boundary condition -B(So) = 1 and those that do not. 

Starting with equation (7), we may set \^ = for a non- rotating source (this is 
a good approximation for slowly rotating sources and weak gravitational fields also). 
Equation (7) then becomes separable, i.e., solutions of it can be written in the form 
B{p,6) = 1 — 2F{p)G{6). The new functions F and G must then satisfy the ordinary 
differential equations 
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(PG „ dG 



0, 



(A.2) 



where (3 is some complex- valued constant. Restricting (3 to be real and requiring that 
/3> — |, the general solutions of equations (A.l) and (A.2) may be written in the form 




"0 



1- 



7:^2 
"0 



(A.3) 



G^O) = C-^i^i^.i (cos^) + C^±Q±^_i (cos^), (A.4) 

where Pl^i^x), Q'^i^x) are the usual associated Legendre functions of the first and second 
kind, respectively, and where c^^^, cj^^, Cg^ and are (dimensionless) constants. Note 
that the solutions (A.3) and (A.4) arc real-valued functions, and so arc P^(a;), Q'^{x) 
(for real fi, v) since they arc defined on the cut (—1, 1) by averaging the relevant limiting 
values of the corresponding complex- valued quantities Pl^{z), Qy{z) [8]. 

We will now find solutions (A.3) admitting the boundary condition F^±(So) = 0. 



1 I r /2Ho(i-f^) 

To do that, we first notice that, since PKyJl — ^) — y "° and since only the 

trivial constant solution is obtained from Pj ^ (y^l — ^) = \J^^-, choosing the -I— sign, 
with Co_,_ = 1 and a suitable choice of Co_,_, the value ^ = corresponds to the unique 
spherically symmetric solution found in ref . [4] . That solution is unique since the function 

_i 1 1 I 5" 

Q 1 ^ {x) differs from P^ {x) only by a numerical factor, and since QlU 1 — ^) again only 

2 2 2 V "0 

yields the trivial constant solution. Besides, since the function Qq{cos6) has singularities 
whenever ^ = or ^ = tt, the corresponding solution cannot be physical (but might be 
considered as a mode solution). Moreover, P_i(cos^) = Po(cos^) gives nothing new, and 
since Q-i{cos9) is undefined, that expression cannot be considered at all. 



Secondly we notice that, since choosing /3 — 2 yields PK^l — = — and 
Pih^l - = y5(^arccos(^l - |^) - which do not fulfil the required 

3 _3 

boundary condition, and since Ql{x)=0 on the cut {Qi^{x) is undefined), the next 



_5 

suitable solution is found by choosing ^ — 6. This is so, since even though Pj ^ (, / l — |^ 

2 V ^0 ' 

1 I T /2Sg(l-4) 

does not fulfil the required boundary condition, Pi^(wl — ^) = 3^ — ° obviously 
does. This solution was found in ref. [5] and corresponds to a pure quadrupole field since 
it involves the Legendre polynomial P_3(cos^) = P2(cos^) (also, Q2{cos6) is singular and 
the corresponding mode solution thus unphysical, and Q-^{cos,6) is undefined). Moreover, 
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given the required boundary condition, this solution is also unique since Ql{x)=Q on the 

2 

_ 5 

cut {Qi '^{x) is undefined). 

2 

Similarly, since the subsequent suitable solution is found by choosing /3 = 20 and 
the +-sign, corresponding to a pure octopole field involving the Legendre polynomial 
P_5(cos6') = P4(cos6'), it would seem that all suitable solutions are given by choosing 
the +-sign and (3 — 2n{2n + 1), n = 0,1,2,.., corresponding to pure even multipole 
fields involving the Legendre polynomials P_2n-i(cos^) = P2n(cos^^). With the required 
boundary condition, these solutions are also unique, since the functions P^"^ ^(^/l — 

2 V "0 

Tn-\- — 

m — 0,1,2,.., do not admit it, and since Qi ^(a;)=0 on the cut (found by using the 

2 

recurrence relation Qi^+^{x) + 2{n+l)x{l-x^y^/^Q'^^^{x) + {u-n){u + n + l)Q'^{x) = 0). 
Also, Qi ^(x) and (5-m-i(cos^) are all undefined. Finally, since Q2n{cos9) is singular 

2 

for all n, the corresponding mode solutions are unphysical. 

That the above extrapolation is correct, can be checked by using the formulae [8] 



r 4± 



,^V^(0)=-2^^-V^sm[|(i.;^)]|^^ (-) 

Since the reciprocal gamma function l/r{x) possesses simple zeros at x = 0, —1, —2, .., 
[8], these expressions vanish for the chosen values (3 = 2n{2n + 1), n = 1,2,3,.., and 
using the +-sign (the special case /3 = was addressed previously). Moreover, no other 
values of /3 will do (since the other possible values f3 = 2n{2n — 1) yield functions of the 
form Qi ^ (x)=0), so we have really found all relevant solutions satisfying the boundary 

2 _ 

condition P^±(So) = (using the — sign rather than the +-sign yields no further zeros). 

So, to summarise, said boundary condition is fulfilled by choosing the functions 
Pi ^(^/l — ^) in equation (A. 3); all other choices fail (or are irrelevant). 

2 V "0 
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